Detecting Quantum Critical Points using Bipartite Fluctuations 
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We show that the concept of bipartite fluctuations T provides a very efficient tool to detect quan- 
tum phase transitions in strongly correlated systems. Using state of the art numerical techniques 
complemented with analytical arguments, we investigate paradigmatic examples for both quantum 
spins and bosons. As compared to the von Neumann entanglement entropy, we observe that T 
allows to find quantum critical points with a much better accuracy in one dimension. We further 
demonstrate that T can be successfully applied to the detection of quantum criticality in higher 
dimensions with no prior knowledge of the universality class of the transition. Promising approaches 
to experimentally access fluctuations are discussed for quantum antiferromagnets and cold gases. 
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Quantum phase transitions [1] occur at zero temper- 
ature and are solely driven by quantum fluctuations. 
Hence it is expected that a quantum phase transition 
should be manifested through the system's entanglement 
properties [2]. Identifying appropriate measures of en- 
tanglement is, however, a non-trivial task. An important 
tool to access and quantify the amount of entanglement 
between two sub-sets A and B of an interacting quantum 
system is the von Neumann entanglement entropy (EE). 
In one dimension (ID), conformal field theory and exact 
calculations have established the logarithmic scaling of 
the von Neumann entropy [3] for critical systems. For 
gapped systems the EE saturates to a constant and thus 
obeys a strict area law (assuming a local Hamiltonian) [4]. 
In fact, EE can help to locate the quantum critical point 
(QCP) in some cases [5]; for more subtile situations (e.^. 
like Kosterlitz-Thouless transitions) it was demonstrated 
recently that the EE failed to locate the QCP of the 
frustrated J1-J2 chain [6]. In higher dimensions, it was 
established that the gapless Heisenberg antiferromagnet 
(AF) on a square lattice obeys a strict are law [7, 8], as 
also expected for a gapped phase. In such a situation, 
it is therefore unlikely that von Neumann EE will be a 
useful and practical tool to detect QCPs. Conversely, 
the valence bond entropy has been shown to be a pow- 
erful quantity to locate QCPs in any dimension, based 
on different scaling regimes, but it is restricted to SU(2)- 
invariant spin systems [6, 9]. 

The aim of this Letter is to promote a general and more 
practical quantity to precisely locate QCPs for a larger 
variety of strongly correlated systems in any dimension 
d. Using the concept of bipartite fluctuations [10-12] 

of particle number or magnetization in many-body 
quantum systems, we focus on systems where such U(l) 
charges O are globally conserved while they locally fluc- 
tuate within each subsystems. We define 

ieA ieA 

where the (globally) conserved quantity O can be the 
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FIG. 1: (color online). For a c/- dimensional system, the fluc- 
tuations T within subsystem A (of linear size L) with respect 
to B provide a precise estimate to locate a QCP at Ac be- 
tween gapless (quasi) ordered and gapped disordered phases 
with disctinct scalings with L. 



particle number n or the magnetization and (•) refers 
to the ground state at T = 0. Oi is defined for a sub- 
system A embedded in a larger one, see Fig. 1. For 
the special case that A is the total system, J^a is just 
the susceptibility (or compressibility, respectively) di- 
vided by temperature. We show for various models, 
such as the spin-^ frustrated Ji — J2 AF in ID, the 
Bose-Hubbard chain at unit filling, 2D coupled Heisen- 
berg ladders, and Bose-condensed hard-core bosons on 
a square lattice, that provides a very efficient tool 
to accurately detect quantum criticality in the frame- 
work of quantum Monte Carlo (QMC) and Density Ma- 
trix Renormalization Group (DMRG) simulations on fi- 
nite size systems [13]. The key feature of the bipartite 
fluctuations is the distinct scaling behavior for gapless 
and gapfull phases in any dimension d [10, 11], as sum- 
marized in Fig. 1: within a sub-system of linear size L, 
exhibits a strict area law for a disordered (gapped) 
ground-state, J-gapped L^~^ ^ whereas for a (quasi) or- 
dered gapless state multiplicative logarithmic corrections 
appear, J-gapiess ^ L^^^lnL, thus allowing to precisely 
locate a QCP between two such regimes. The bipartite 
fluctuations give an alternative view of the correlation 
functions since they are dominated by short-range fluc- 
tuations [11]. Experimentally, the concept of fluctuations 
has a very strong potential [12]. 

One dimensional systems — We now address ID mod- 
els, governed by Kosterlitz-Thouless (KT) type quantum 
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FIG. 2: (color online). Luttinger parameter K of the Ji — J2 
chain extracted via (3) vs. A = J2/J1. Shown is DMRG data 
for L = 100 (red squares) and 200 (black dots) for PBC. 

phase transitions usually difficult to precisely locate nu- 
merically. The first model we study is the frustrated 
spin-^ J1-J2 chain, governed by the Hamiltonian 

H(A) = ^(S,-Si+i + ASi-Si+2) , (2) 

i 

where J2/J1 = A > 0. For A < Ac, this system has 
power-law critical correlations. At Ac — 0.2412, a KT 
transition into a dimerized phase occurs [14, 15]. As men- 
tioned above, the estimated value for the QCP using EE 
is not very precise compared to the established meth- 
ods [14] because the prefactor of the leading term in the 
EE (i.e., the central charge c) is more or less insensi- 
tive to a change of A close to the QCP [6]. Instead, we 
detect the transition by observing the behavior of un- 
der variation of the control parameter A which triggers 
the quantum phase transition. The low-energy theory 
describing such a quasi-ordered state is the Tomonaga- 
Luttinger liquid [15], for which [11, 12] 

K 

T{L) = ^InL + cst, (3) 

where i^^ = 1/2 is the Luttinger liquid parameter of this 
SU(2) point. However, marginally irrelevant operators 
lead to sizeable logarithmic corrections for K [16], when 
computed on finite size systems. Interestingly, such cor- 
rections vanish precisely at Ac where K quickly reaches 
its asymptotic value of 1/2. Thus we have a systematic 
method at hand to detect this phase transition. In Fig. 2 
we have plotted the Luttinger parameter K extracted 
from finite size DMRG calculations of Eq. (3) versus A. 
For PBC and L = 100, 150, 200, and 250, and after 
performing finite-size scaling, we obtain Ac = 0.2412(3) 
which agrees very well with the best estimates [14] . While 
there are a few other techniques available to find the QCP 
of the J1-J2 chain [6, 14, 17-19], our approach stands out 
through efficiency and simplicity. 

Another interesting model is the Bose-Hubbard chain: 

H = -tJ2blb, + ^^ni{ni-l)-J2t^ni , (4) 



. L = 256 (OBC) 




£ 1.8 






g 0.294 0.298 0.302 



0.3 0.4 4- IT J 0.5 



FIG. 3: (color online). Luttinger parameter K of Bose- 
Hubbard chain extracted via (3) vs. A = t/[/ for L — 256 
(unit filling) and OBC. We restricted the local boson occupa- 
tion number to 4 [21, 22]. Inset: zoom close to the transition. 

where t is the hopping amplitude, U the on-site repul- 
sion, and \i the chemical potential. Away from half fill- 
ing, we expect a superfluid-Mott transition triggered by 
A = tjU . The superfluid phase is a Luttinger liquid [20] 
with Luttinger parameter K > \. For unit filling, the 
QPT from a superfluid to a Mott insulator is of KT type 
(like in the J\ — J2 chain discussed above). The com- 
plete {ji^t/U) phase diagram was carefully investigated 
within DMRG in Refs. 21-23. Here we revisit the prob- 
lem (restricted to unit filling) and show that we locate 
the transition with a better accuracy by virtue of the fluc- 
tuations. In the superfluid phase, the Green's function 
G(r) = {h^ho) (X r~^^'^^ decays as a power-law. From 
Luttinger liquid theory we know that the transition oc- 
curs for Kc = 2, see Ref. 24. In Refs. 21, 22 the Luttinger 
parameter K was extracted directly from G(r), thus giv- 
ing an estimate of the critical point Ac = 0.297±0.01 [22]. 
The major advantages of our approach is that (i) we have 
a finite size formula for the fluctuations (i.e., applicabil- 
ity of conformal mappings) contrary to G(r), and (ii) 
the computational cost of using DMRG (see App. C of 
Ref. 12) is much lower as compared to the Green's func- 
tion at large distances. We extract K from for OBC 
with L = 64, 128, and 256 (the latter is shown in Fig. 3). 
By performing finite size scaling we obtain a much more 
precise estimate Ac = 0.2989(2), as compared to previous 
works [13]. 

Two dimensions — Let us now move to 2D with a sys- 
tem of coupled spin-^ AF ladders, depicted in the inset 
of Fig. 4 (a), and governed by the Hamiltonian 

U=Y, ^i-^j + E ^S,-S,. (5) 

ladd. inter— ladd. 

This model [25, 26] displays a gapped valence bond solid 
(VBS) phase for small inter-ladder coupling A < Ac with 
Ac = 0.31407(5) [25], and a gapless Neel ordered phase 
for A > Ac- Here we investigate the T = fluctuations 
of the total magnetization in a region A of size x x y 




FIG. 4: (color online). Quantum Monte Carlo results for T = fluctuations T of the total magnetization in a region A for 2D 
coupled spin- 1 ladders [Eq. (5)], depicted in the inset of (a). Left (a): T jL increases logarithmically with L in the Neel regime 
(black squares A = 1) whereas it saturates to a constant in the valence bond state (green circles A = 0.1). Right (b): J^/L, 
plotted vs. A for various system sizes, displays a crossing point at Ac. Insets: (%) crossing of the stiffness x L at Ac for the 
same sizes; 1/ L convergence of the crossing point for T (red squares) and ps (black circles) to the critical value (horizontal 
black line) Ac = 0.31407 [25]. 



embedded in a periodic square lattice L x L. We choose 
a sub-system A with x = L/2 and y = L which contains 
an even number of sites. QMC results for the T = [27] 
expectation of J^{L/2) are shown in Fig. 4, with square 
lattices size up to L x L = 10^, for the isotropic square 
lattice A = 1 (Neel) and for weakly coupled ladders with 
A = 0.1 (VBS). In contrast with the entanglement (or 
Renyi) entropy which displays a strict area law in the 
Neel phase [7, 8] (and presumably also in the VBS phase), 
the fluctuations follow a rather different scaling [8]: 
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Therefore, T ji plotted for different sizes will display a 
crossing point at Ac, as we indeed observe in the panel 
(b) of Fig. 4 where the curves T{L/2)/L are plotted for 
various system sizes. The spin stiffness p^, also known to 
be a useful quantity to locate a QCP, is shown in the right 
inset (ii) of Fig. 4 (b) where one sees a similar crossing for 
ps X with z = 1 and d = 2. As usual for such a 

technique, a drift of the crossing point is observed with L, 
as visible in the left inset (i) of Fig. 4 (b). Already known 
for a few other models [26, 28], the crossing points ob- 
tained from the stiffness converge very rapidly with 1/L 
to the bulk value Ac, whereas we found a slower conver- 
gence for the estimates obtained from T/L. Despite such 
effect (which may not be generic but model dependent), 
we demonstrate here with this simple example that is 
a very useful quantity to locate a QCP between ordered 
and disordered phases for d > 1. 

One can get even more insight from the behavior of 
the coefficients a and P in Eq. (6) as a function of the 
inter-ladder coupling A (see Fig. 5). The prefactor a 
of the leading term ~ LlnL in the Neel phase van- 
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FIG. 5: Prefactors a and /S from Eq. (6) for coupled Heisen- 
berg ladders, extracted from QMC data of Fig. 4, and plotted 
against A. (a) The critical point is shown by a red circle, and 
the green curve is the power- law fit indicated on the plot, (b) 
The vertical dashed line signals the critical coupling Ac and 
the crossing point (red circle) is at /Sc — 0.0835. 



ishes at the critical point a ^ {X — Ac)^, with x 0.7 
and Ac = 0.315(1), in good agreement with the value 
0.31407(5) [25]. The area law term displayed in 
Fig. 5 (b), although certainly non- universal, exhibits a 
very interesting A-shape and passes through a maximum 
/3c — 0.0835 at the critical coupling Ac- 

It is important to emphasize that, contrary to the stiff- 
ness, a prior knowledge of any critical exponent, such as 
the dynamical exponent is not necessary to precisely 
locate the QCP. Note also that we expect the valence 
bond entropy [9] to display similar crossing properties for 
such a SU(2) symmetric Hamiltonian Eq. (5). In order 
to illustrate further the general character of this method, 
we focus now on a non-SU(2) model: hard-core bosons 
on the square lattice. Governed by the Hamiltonian 
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where b are hard-core bosons operators, t the hopping 
integral and fi the chemical potential, hard-core bosons 
on the square lattice [29] exhibit a particle-hole symmet- 
ric phase diagram at T = with a Bose-condensed su- 
perfluid state for < 4, and trivial Mott insulating 

phases for \ fi/t\ > 4, the transition between them being in 
the universality class of the diluted Bose gas with z = 2. 
The Bose-condensed (U(l)-broken-symmetry) state is ex- 
pected to display for (fluctuations of the particle num- 
ber) a similar scaling as the one observed for SU(2)- 
broken Neel-ordered spins, whereas for the trivial Mott 
insulators we simply have J-Mott = 0. In Fig. 6 (a) T = 
QMC results obtained for are shown for 4 representa- 
tive values of the chemical potential. The prefactor a of 
the L In L term is plotted versus the chemical potential ja 
in the right panel of Fig. 6 where we observe a very inter- 
esting dome-like shape in the superfluid regime. One can 
use an interesting analogy with quasi-one dimensional 
systems where the Josephson type interchain tunneling 
term will lock the superfluid phase difference between all 
chains. The low-energy (quasi-ordered) superfluid phase 
is described in terms of a single macroscopic ID gapless 
mode. For a number of chains N = L then we pre- 
dict = {KL/tt'^) InL. The logarithmic scaling of is 
controlled by the Luttinger parameter K of the effective 
theory. In the hydrodynamic description of a Luttinger 
liquid K = TrV^A^Tsf , where is the compressibility and 
Tgf is the stiffness. This gives a = V^Tgf/Tr. A similar 
quantum-hydrodynamic theory for interacting bosons is 
obtained in two dimensions using the Gross-Pitaevskii 
approach. Comparing the prefactor a with ^/hiTs^ (ob- 
tained in the same QMC simulation), as shown in Fig. 6 
(b), gives a very good agreement. We find the following 
result for the entire superfluid regime: a(/i) = V^Tgf/p 
with a coefficient p 3.2(1). Scaling relations close to 
a QCP at Ac predict Tgf - ^ ^ ^z-d^ ^j^^g 

leading to a ~ ~ ^c)^, which can be compared to 
Fig. 5 (a) where the exponent x 0.7 is very close to 
v = 0.709(6) of the 3D Heisenberg universality class [26]. 

Conclusion — The concept of bipartite fluctuations of 
a (strongly correlated) quantum system has been shown 
for various paradigmatic condensed matter models to be 
an efficient, accurate, and rather general tool to detect 
quantum critical points using state of the art numerical 
techniques. In contrast to the von Neumann entropy, the 
fluctuations can be successfully used even in two spatial 
dimensions to find the critical point. Promising paths 
to directly measure the fluctuations have been proposed 
recently [12]; particularly interesting proposals are quan- 
tum magnets in an external magnetic field with Meissner 
screens (covering region B) as well as direct measurement 
of using single atom microscopes [30] . A next step will 
be to test the usefulness of this tool for unconventional 
quantum criticality [31]. 
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FIG. 6: (color online). QMC results for the fluctuations T 
of the particle number for hard-core bosons on square lattices 
of size L X L, the subsystem A being half the lattice (see in- 
set), (a) J^/L plotted against InL for 4 representative fillings 
indicated on the plot, and compared to the SU(2) Heisen- 
berg model, (b) Prefactor a, obtained from fits to the form 
Eq. (6), shown over the entire regime versus /x/t together with 
V^^i/p from QMC, with p = 3.2(1). 
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Supplementary Material for 
"Detecting Quantum Critical Points using Bipartite Fluctuations" 



Bipartite fluctuations [1] bring a very strong concept in 
the study of quantum criticality. Our approach to locate 
Quantum Critical Points stands out through the combi- 
nation of the attributes generality, simplicity, efficiency, 
and accuracy. While other existing methods might be 
comparable with regard to one of these attributes, none 
of them is comparable with regard to all of these at- 
tributes as we shall show in the following. While it would 
be an impossible task to prove this statement for any 
known quantum critical point within any numerical ap- 
proach, we rather focus here on a representative exam- 
ple and leave it to the reader to convince himself/herself 
about the usefulness of our method. We show below 
for the Kosterlitz-Thouless transition between superfluid 
and Mott insulator phases of the Bose-Hubbard chain 
at unit filling (numerical hard to study due to finite size 
logarithmic corrections) that the comparison between our 
approach and several existing ones is unambiguously in 
favor of our method. 

GENERALITY 

Our method is applicable to spins and interacting 
bosons (both demonstrated in the main text) but also 
to interacting fermions. It works equally well in any di- 
mension (demonstrated for d = 1,2 in the main text). 
Apart from the U(l) symmetry of spin or charge (which 
is of course essential for our approach), presence or ab- 
sence of any symmetry such as e.g. SU(2) does influence 
the precision of our method. Moreover, no prior knowl- 
edge about the order parameter or any critical exponent 
is required to locate the quantum critical point. Finally, 
our approach even enables us to locate quantum critical 
points of Kosterlitz-Thouless type (as demontrated for 
two different models in d = 1), notoriously cumbersome 
due to finite size logarithmic corrections. 

SIMPLICITY 

The fiuctuations can be easily computed for a given 
subsystem, numerically and often even analytically. For 
numerical approaches where the reduced density ma- 
trix of a subsystem is computed the bipartite fiuctu- 
ations are a side-product (see App. C of Ref. 1). Ex- 
amples are the density matrix renormalization group 
(DMRG) and its descendants. For other numerical meth- 
ods, J^A requires to compute diagonal correlation func- 
tions G{i,j) = (S^Sj) for magnets or {niUj) for itin- 
erant systems. Simple crossing technique, i.e., plotting 

I j^d-i ^ different system sizes L reveals the 
quantum critical point. 



EFFICIENCY 

Our approach is very efficient for the DMRG technique 
where the reduced density matrix can be directly com- 
puted. For Quantum Monte Carlo techniques, while J^a 
is not conserved during the propagation along the Trot- 
ter direction, one can easily keep track of its evolution 
without computing all correlators at each imaginary time 
slice. As a result, the computational cost for estimating 
F is the same as of getting either the stiffness or a given 
structure factor. If the type of order is unknown, the full 
set of correlations will be required, which is, computa- 
tionally speaking, much more expensive. In such a case, 
computing would be far much easier. 

Alternatively, one could compute the gap for various 
system lengths, but then one needs to perform a finite 
size scaling analysis which involves a prior knowledge of 
the universality class of the transition. Moreover, an ef- 
ficient estimate of a tiny gap with QMC or DMRG (also 
discussed below) is computationally more costly than the 
other observables discussed above. 

If one knows which kind of order the system will 
achieve, the natural computation would be the associ- 
ated structure factor, with the complication that finite 
size scaling at the transition may be tricky to analyze. 
Therefore one usually prefers to rely on crossing tech- 
niques. For instance the order parameter squared (struc- 
ture factor) S X L^^/^ will displays a crossing at a QCP, 
but one needs to know the ratio 13/v oi two important 
critical exponents. The stiffness ps is usually better since 
it only requires to know one exponent: the dynamical 
exponent z. The good thing with the crossing of the 
stiffness is that subleading finite size corrections are usu- 
ally very small (see for instance Ref. 2) and therefore the 
location of Ac is very good. Nevertheless, in some case 
the dynamical exponent is an important unknown quan- 
tity of the problem and the crossing of ps x L'^-^-z 
be useless (here d is the spatial dimension). As a con- 
sequence, the J^-crossing technique, which only relies on 
the result that J^{x) ~ x^~^hix, provides a very nice 
and alternative way without any assumption regarding 
the critical point universality class. 

In the following we consider the one-dimensional Bose- 
Hubbard model and use the DMRG method. First, we 
compute the fiuctuations for a given parameter setting 
and, second, the Green's function G(r) = {h\h^) for the 
same setting. We eventually compare the used CPU 
times. 

CPU time depends on many internal {i.e., method- 
specific) parameters of the DMRG method. We tried to 
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FIG. 7: Normalized CPU time using a single core vs. kept 
DMRG states for Bose-Hubbard chain, L = 128, OBC. For 
about 300 kept DMRG states convergency is guaranteed im- 
plying a normalized CPU time > 2. 



use a standard setting which is representative for most 
DMRG implementations. We always computed the fluc- 
tuations and the correlation functions on the same com- 
puter using a single core (in order to rule out the influ- 
ence of the parallelized code). We also kept the same 
parameter setting for all computations. In addition, we 
varied the number of kept DMRG states (which can be 
converted into a measure for convergency) and for dif- 
ferent interaction strengths t/U. Since computations for 
different values of t/U have been performed on different 
computers, we normalized the CPU time. 



normalized CPU time 



CPU time for correlations 
CPU time for fluctuations 



For 300 kept DMRG states we found the "discarded en- 
tropy sum" to be smaller than 10~^^ which guarantees 
convergency within the DMRG method (we always per- 
formed 15 sweeps). For all considered values t/U (see 
Fig. 7) the normalized CPU time is between 2.1 and 2.3, 
ie., it takes more than the double time to compute cor- 
relations compared to fluctuations. In this example, we 
chose the system length to be L = 128 sites and we used 
OBC. In the main text, we also computed L = 256 where 
the normalized CPU time is even higher. Below we fur- 
ther demonstrate that T is the best quantity to easily 
locate the critical point of the ID Bose- Hubbard model, 
as presented in the main text. 

Before this, let us briefly comment on investing the gap 
size as an alternative method to locate the quantum criti- 
cal point. For most numerical methods, the performance 
to compute the gap size between ground state energy 
and excitation spectrum is rather bad. In the DMRG 
method, for instance, the normalized CPU time would 
easily exceed a value of 4 to 5. 



ACCURACY 

For the crossing technique, i.e., plotting T/L^~^ vs. 
A for different system sizes L, one only needs to know 
the value of the fluctuations when subsystem A is e.g. 
half the total system length. No fitting process is in- 
volved (in contrast to correlation functions (see below)). 
The crossing technique works extremely well as demon- 
strated in our paper for both two-dimensional models. 
In one dimension, it works as well, but here in order to 
beat the existing approaches (such as level-spectroscopy 
of the full energy spectrum) we used a somewhat more 
sophisticated way and extracted the Luttinger parameter 
K from the fluctuations. This involves a fitting process 
but due to the availability of conformal mappings we have 
a finite size formula at hand which simplifies the fitting 
process drastically and makes it much less dependent of 
the fitting window. In contrast, extracting K from the 
correlation function appears to be a much more delicate 
procedure. For the Green's function G{r) we still have a 
finite size formula [3] , 



^y sin(7rr / L) sin(7r / L) 



2Ll sin(7r(r + 1)/(2L)) sin(7r(r - 1)/(2L)) , 



but the extracted value of K is very sensitive to the fit- 
ting window which is used. This is exemplified in the left 
panel of Fig. 8 below, where DMRG results are shown 
for an open Bose-Hubbard chain of L = 128 site, in the 
superfluid regime (n = 1, t/U = 0.3), but close to the 
Mott insulating phase. Indeed, the extracted value of 
the Luttinger liquid parameter K strongly depends on 
the fitting process, depending on the type and the size 
of the retained fitting window. Consequently the esti- 
mated K = 1.93(11) displays a quite large uncertainty, 
as already discussed in T. D. Kiihner, S. R. White, and 
H. Monien (Ref. 10). In contrast, the bipartite fluctua- 
tion provide a much more accurate tool, as shown in the 
right panel of Fig. 8 where fitting using the finite size 
formula [4] 



K 
2^ 



In ( — sin(7rr/L) ) + const.. 



leads to much more precise estimate oi K = 2.03(2), 
much less sensitive to the fitting process. Given the fact 
that the numerical computation of is twice less numer- 
ically expensive than G(r), this specific example clearly 
tells us that our method is superior to existing ones. We 
further demonstrate this fact in the following where we 
present a quantitative comparison with several previous 
studies, performed over the last 20 years. 
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FIG. 8: DMRG results for the Bose-Hubbard chain at t/U — 0.3 for L — 128 (OBC) is shown. Retaining the same number of 
DMRG states (300), the required CPU time to compute the Gree's function G{r) (Left) was twice larger that for the bipartite 
fluctuations T{r) (Right). While extracting K from the correlations is very sensitive to the fitting window/range (yielding 
K — 1.93(11)), it is much less sensitive when extracting K from the fluctuations T (yielding K — 2.03(2)). 



CRITICAL POINT ESTIMATES OF THE SF-MI TRANSITION IN THE ID BOSE-HUBBARD MODEL 

AT UNIT FILLING 



Year 


Reference 


Technique 


Observable 


Estimate 


1991 


Krauth [5] 


(approximate) Bethe Ansatz 




1/(2^3) - 0.2887 


1992 


Batrouni et al. [6] 


QMC 


Superfluid stiffness 


0.2100(100) 


1994 


Elesin et al [7] 


Exact Diagonalization 


Gap 


0.2750(50) 


1996 


Kashurnikov et al. [8] 


QMC 


Gap 


0.3000(50) 


1999 


Elstner et al. [9] 


Strong coupling 


Gap 


0.2600(100) 


2000 


Kiihner et al. [10] 


DMRG 


Correlation function 


0.2970(100) 


2008 


Zakrzewski et al. [11] 


Time Evolving Block Decimation 


Correlation function 


0.2975(5) 


2008 


Laiichli et al. [12] 


DMRG 


von Neuman entropy 


0.2980(50) 


2008 


Roux et al. [13] 


DMRG 


Gap 


0.3030(90) 


2011 


Ejima et al. [14] 


DMRG 


Correlation function 


0.3050(10) 


2011 


Danshita et al. [15] 


Time Evolving Block Decimation 


Excitation spectrum 


0.3190(10) 


2011 


This work 


DMRG 


Bipartite Fluctuations 


0.2989(2) 
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FIG. 9: Various critical point estimates of the SF-MI transition in the ID Bose-Hubbard model at unit filling. The plotted 
values in function of time refers to the above table and to the reference list below. 
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